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ABSTRACT 

We combine the globular cluster data for fifteen Brightest Cluster Galaxies 
and use this material to trace the mass-metallicity relations (MMR) in their glob- 
ular cluster systems (GCSs). This work extends previous studies which correlate 
the properties of the MMR with those of the host galaxy. Our combined data sets 
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show a mean trend for the metal-poor (MP) subpopulation which corresponds to 
a scahng of heavy-element abundance with cluster mass Z ~ ^-soio.os^ trend 
is seen for the metal-rich (MR) subpopulation which has a scaling relation that 
is consistent with zero. We also find that the scaling exponent is independent of 
the GCS specific frequency and host galaxy luminosity, except perhaps for dwarf 
galaxies. 

We present new photometry in (g',i') obtained with Gemini/GMOS for the 
globular cluster populations around the southern giant ellipticals NGC 5193 and 
IC 4329. Both galaxies have rich cluster populations which show up as normal, 
bimodal sequences in the colour-magnitude diagram. 

We test the observed MMRs and argue that they are statistically real, and 
not an artifact caused by the method we used. We also argue against asym- 
metric contamination causing the observed MMR as our mean results are no 
different from other contamination-free studies. Finally, we compare our method 
to the standard bimodal fitting method (KMM or RMIX) and find our results 
are consistent. 

Interpretation of these results is consistent with recent models for globular 
cluster formation in which the MMR is determined by GC self-enrichment during 
their brief formation period. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: individual (IC 
4329, NGC 5193) 



1. Introduction and Background 



Colour bimodality in globular cluster systems (GCSs) has been observed for nearly all 



types of galaxies, i ncluding smal 



Harris et al 



2006 



2006 



Mieske et al.l 



Cantiello et al 



12007 



and large galaxies, and elliptical and spiral galaxies (e.g 



2006 



Wehner et al 



Peng et al.l l2006al : ISpitler et all 12006 



Strader et al. 



20081 ) and is now thought of as a "universal" char- 



acteristic that reflects their formation history. 

This red/blue bimodality corresponds to a bimodality in metallicity. Concern has 
been raised that the split occu rs because of a non-linear colour-metallicity relation (e.g.. 



Worthevlll994 



Yoon et al. 



IR colours (IKundu fc Zepi 



2006); however, more dir ect metallicity measurements, including 



20071 ) and spectroscopy (IBeasley et al.ll2006l : IStrader et al.l 12007 



Beasley et al. 20081 ) . argue against this. In addition, t he GCS for the M ilky Way is found 



to have a linear (^i?-/j-metallicity relation (Figure 7 in lHarris et al.ll2006l ). M31 also shows 
bimodality in (U-V)o, (U-R)q and (V-K)q colours and metallicity, through spectroscopic and 
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photometric studies presented in iBarmby et al.l (|2000| ). Bimodality for M31 is not seen in 
othe r colours, inclu ding (V-I)o - a colour often used previously with the testing of bimodal- 
ity. iBarmby et al.l suggest, however, that photometric and reddening estimate errors blur 
away the expected bimodality distribution in these colours (that are actually less sensitive 
to metallicity). 

Another concern is that age and metallicity are degenerate. One attempt to understand 
this issue was undertaken by iPuzia et al.l (120051 ) who used Lick indices to look at the GCSs 
of seven early-type galaxies. They found that there are three age groups into which GCs 
fall: two-thirds have an age >10 Gyr, one-third are between 5 and 10 Gyrs old, and a 
small fraction have an age <5 Gyrs. The issue of age- metallicity degeneracy remains to be 
resolved, but the spectroscopic data provide clear evidence that colour bimodality transforms 
to metallicity bimodality. If age alone could account for the GCS bimodality, a difference 
of A{B — I) ^ 0.3 would correspond to age differences of 4.5 Gyr and 7.5 Gyr for constant 
metallicities of f Fe/H]=-2.25 an d -1.5, respectively, with the redder clusters being the older 



subpopulation (IMarastonl |2005| ) . 



Remaining challenges to unravel are the presence of the two distinct subpopulations 
and what these may r aean for galaxy f ormat ion theories. The ideas can be grouped into 
three main categories: lAshman &: Zepj (119921 ) suggest gas-rich galaxy merging, where each 
progenitor galaxy has its own blue, metal-poor (MP) GCS and then upon merging the two 
galaxies create a re d, metal-rich (MR) subpopulation. In multi-phase dissipational collapse 
(iForbes et al.lll997l ) the GC forma tion process is halted, for example , via cosmic reionization 
after a galaxy forms its MP GCs (iBeaslev et al.ll2002l: ISantosl 120031 ) and then resumes some 
time later to form MR GCs. Finally, ICote et al.l (119981 ) suggest that an initiall y large galaxy 
has its own MR GCs and then accretes MP GCs from satellite galaxies. iHarrid (120031 ) 
comments that these processes are not mutually exclusive, and are likely to act in a number 
of combinations for different circumstances. 

In addition to the bimodality, a correlation between luminosity and mean colour - that 
is, a mass-met allicity relation (MMR) - is seen for the blue subpopulation iii at least some 



mass ive galaxies (e.g.. lHarris et al.ll2006l : iMieske et al.ll2006l : ISpitler et al.ll2006l : IStrader et al. 
20061 ). This feature is also referred to as the "blue tilt" . Drawing material from the previously 



published literature, we list values in Table [T] of the measured MMR slope along the blue GC 
sequence; this value, p, is the exponent in the scaling of heavy-element abundance Z with 
GC mass M, Z^M^. We will discuss the derivation of these slopes and their comparisons 
in Sections 3 and 4. Columns in Table 1 contain 1) the galaxy sample from which the blue- 
MMR value was obtained; 2) the value of p for the blue-MMR with errors where available; 
3) the brightest GC luminosity values that were used to make the fit, and 4) the reference 
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to the data source paper. 



No obvious MMR al o ng; the red sequence has been obs erved (Harris et aLll2006l : iMieske et al, 
20061 : ISpitler et allbood : IWehner et allboosh . However, IWehner et all J2008h detected the 



red sequence continuing to brighter magn i tudes towards the ultra-compact dwarf (UCD) 
galaxy range. In addition, IStrader et al.l (120061 ) interestingly find no blue-MMR in M49 
(NGC 4472). If confirmed, this result will provide an additiona l comp lexity to be explained 
by galaxy and globular cluster formation models. iMieske et al.l (120061 ) suggest that environ- 
mental causes may produce this result because M49 is at the centre of its own sub-cluster 
within the Virgo Cluster. This may be an example, they claim, of the differences amongst 
galaxies with different accretion/merger histories if MMRs are created through GC accretion 
from galaxies with different masses. Testing this hypothesis, they find that in general GC 
accretion is unlikely to account for the entire slope of observed MMRs, although it may 
contribute to or, in the case of M49, negate the trend. 

A few other galaxies with smaller GC populations have also been searched for an MMR. 
Both red and blue s lopes for the GC sequences in NGC 1533 are consistent with zero 
(IDeGraaff et al.ll2007l ). This SBO galaxy is neither in a cluster nor is it a bright elliptical, 
which the authors imply would be the necessary conditions to produce the self-enrichment 
that they suggest causes the blue-MMR. The Milky Way does not show any evidence for a 
blue-MMR, b ut this may be due to small-number statistics - as there are only ^15 GCs 
(iHarrid Il996l ) - rather than the fact that it is a spiral galaxy. ISpitler et al.l (12006! ) show 



that blue-MMRs are not limited to ellip ticals as they see a blue-MMR in the Sombrero 
Gala xy (NGC 4594 ) iSpitler et al.l l2006l classify this galaxy although they note 

that iRhode fc Zepj (120041 ) find it to be an SO. As with ellipticals, no MMR is seen for 
NGC 4594's red subpopulation. A key point in assessing the presence or abs ence of a de 



tectable slope has proven to be the total cluster population in the galaxy (e.g., IMieske et al. 



20061 ) : if the MMR slope is most noticeable for the high-luminosity end of the GC sequence 
(M > IO^Lq), it will be unambiguously detectable only in galaxies with thousands of clusters 
where the high— L end is well populated. For most galaxies, and certainly for most spirals 
and dwarfs, this upper end is thinly populated and the existence of the MMR is likely to be 
not decidable. 



Spitler et al.l (120061 ) claim that the blue-MMR is not due to the addition of a group of 
another kind of objects at brighter magnitudes (around My ^ -H); but that the relation 
continues down to fainter magnitudes - at least to the GC lurn i nosity function turnover 
magnitude of A^^'^ = -7.6 ± 0.06. By contrast, iHarris et al.l (120061 ) in their study of 



eight giant ellipticals find that the MMR slope is clearly nonzero only for the brightest blue 
GCs and that the blue sequence is essentially vertical (i.e., no MMR) for L < 6 x lO^L©. 
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Mieske et al.l ( 120061 ) use a sample of several dozen Virgo galaxies, binned into four groups 
by galaxy luminosity, to find that the lowest-luminosity members (the dwarfs) show little or 
no MMR; however, these are also the galaxies for which the high-luminosity end of the GC 
sequence is least well-populated. 

Though this is still a new feature of GCS systematics, several model interpretations 
for the MMR have already been raised. The most likely ones - both observationally and 
statistically - appear to have organized themselves around the early enrichment history of 
GCs during their formation, either from their own internal self-enrichment or from their lo- 
cal surrou nding host ga s cloud s. Working within the hierarchical-merging picture of galaxy 
formation, [Harris et al.l (120061 ) suggest that the blue, MP GCs form first while they are still 
within dwarf-sized and metal-poor host clouds. The more massive ones can enrich further 
during formation, thus giving rise to an MMR. By contrast, the red-sequence GCs then form 
later and are slightly more metal-enriched and red because they form in a de e per p oten- 
tial well when the galaxy formation process has been completed. IStrader et al.l (120061 ) and 
Mieske et al.l (120061 ) favor self-enrichment as th e cause of the blu e MMR . Q uantitative mod- 
els bu ilt on these ideas have been developed by lStrader &: SmithI (120081 ) and lBailin &: Harris 
J2009h . 



Although rare. Brightest Cluster Galaxies (BCGs) have uniquely large samples of GCSs 
which make them the best environments for studying MMRs. In this paper we present 
newly reduced GC photometric data for two BCGs, NGC 5193 and IC 4329. We follow with 
a new ana lysis of the HST d ata fo r four more BCGs a nd further combine them with the 
data from iHarris et al.l (120061 ) and IWehner et al.l (120081 ). We view our discussion in total 
as one more step towards understanding this intriguing new correlation, which promises to 
uncover additional clues towards GC formation and enrichment history. Other papers by 
the authors, both on the observational and theoretical side, are in progress which will help 
develop a more complete synthesis. 

Section [2] details the observations and data reduction of these two galaxies, and is where 
we also obtain the colour-magnitude diagrams (CMDs) and analyse the MMR. We compare 
our work with other similar studies and discuss the results in Section [31 we summarize our 
work in Section HI 



2. New Observations and Data Reduction 



We first discuss our new observations for the giant galaxies NGC 5193 and IC 4329. 
Some uncertainty exists about the status of NGC 5193 and the cluster of galaxies of which 
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it is a part. At times it has been associated with Abell 3560 (e.g. lAbell et al.l 119891 ). How- 
ever, its recession velocity and its projected position indicate that it is neither found at the 
centre o f A3560 nor is it asso ciated with that cluster. The redshift of A3560 is 14470 ± 123 
km s-i Jwillrnere^aDll999h wher eas NGC 5193 and its companion (NGC 5193A) are at 
3600 km s-i JVettolani et al.lll99oh . and thus NGC 5193 and NGC 51 93A appear as fore 



ground obiects. NG C 5193 and NGC 5193A are located in Abell 3565 (IWillmer et al. 



1999 



BardeUi et al.ll2002[ ). Abell 3565, which has a veloc ity of 3586 ± 45 km s"^ flWillmer et al. 
1999|) and a richness class of 1 JAbell et all Il989|) is a cluster centred on IC 4296, the 
brightest galaxy in the cluster (jVettolani et al. [l990 : Willmer et al. 1999). 



IC 4329 is associated with Abell 3574 which is also kno wn as Shapley 1346-30, K27, the 
IC 4329 group and Centaurus North (Phillipps et al.lll993l and references therein). IC 4329 
also has a richness class of 1 (lAbell et al.l Il989h and a velocity of 4808 ± 21 km s~^ (the 
CMB reference frame; NASA Extragalacitic Database). Abell 3574 is in turn a member of 
the Hydra- Centaurus Supercluster. 

These two galaxies were observed with the Gemini-South CMOS camera in {g', i') on 
the nights of 2006 February 8 and March 23. Images have a 5.5' x 5.5' field of view, an image 
scale of 0.146" per pixel, and come from the data set GS-2006A-Q-24. The two fields can 
be seen in Figure 111 Th ese images were taken during the same run as NGC 3311, published 



m 



Wehner et al.l (120081 ) whose data we re-analyse here to compare with the data from NGC 



5193 and IC 4329. Raw exposure times were 13 x 500 sec in g' and 12 x 500 sec in i', similar 
to the data for NGC 3311. We used the GEMINI/GMOS package within IRAF to preprocess 
all the images for NGC 5193, IC 4329, and used standard-star exposures taken during the 
run to calibrate the data. Adopted values of redshift and luminosity are taken from NED. 
The foreground extinction values were evaluated by using the approximate 1/ A relation and 
the extinctions in the UBVRI system as listed in NED. 

NGC 5193 and IC 4329 are at distances of 56.6 ± 3.8 Mpc and 72.4 ± 4.6 Mpc (calculated 
from their values of cz, in the CMB reference frame, as in NED along with Ho = 70 km s^^). 
For a typical GC half-light radius of 3 pc, a GC at a distance of 60 Mpc would subtend 
an angle of 0.021" which is ~20 times smaller than the ground seeing FWH M, and so GCs 
appear entirely star-hke. Thus, we can discount the possibility suggested by iKundul (120081 ) 
that the blue-MMR is only seen because of aperture-size biasses in high-resolution data. 

To prepare the images for photometry, we median-subtracted the smoothed isophotal 
light, allowing the faint star-like objects to be seen more easily in to near the galaxy centres. 
For the photometry, we used DAOPHOT (th e standalone version 4) and ALLSTAR in the 
same way as described in lWehner et al.l (120081 ). Two of the nine parameters that ALLSTAR 
returned - the magnitude of the object and the associated error - proved useful to identify 
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false "hits" which were clearly associated with the higher background noise near the brightest 
regions. These were removed from all later analysis. Examples of these false identifications 
are shown in the dashed boxes in Figure [2J 

We used DAOMATCH to match the detected objects across the g'- and i'-band images, 
with DAOMASTER refining the estimates. Only x- and y-shifts were allowed, leaving co- 
ordinate shifts accurate to within ±0.1 pixels. Finally, we also used IRAF's GEOMAP as 
a consistency check. The final photometry files of matched objects measured in both {g' , i') 
comprise 1932 objects for NGC 5193 and 3558 for IC 4329. 



2.1. Calibration 



The photometric zeropoints were set from short exposures of standard stars from lLandolt 



( 1992 ). with {g',i') magnitudes from the catalog by E. Wehner. The CMOS calibration 



images of these were obtained on the nights of 2006 March 25, 26 and 27, the same nights 
as the science images. Because only a handful of standard-star exposures were made during 
the two nights, we view these calibrations as strictly preliminary and note that we intend 
to improve them in a future paper through a specially designed photometric observing run. 
Fortunately, as will be seen in the discussion below, our analysis of the MMR and the bimodal 
sequences depend only on the relative colours, not the absolute scale. The peaks of blue and 
red subpopulations for many galaxies are virtually constant at a value of {g' — i')o ~ 0.8 



and (q' — i')o ~ 1.05 (e.g., for M87 with extensive calibration to the SDSS system; [Harris 



20091 ). Using this information, we estimate our colour terms to be accurate for NGC 5193, 



and offset by {g' - i')o ~ 0.2 for IC 4329 and NGC 3311. 



2.2. Colour-Magnitude Diagrams 

The colour-magnitude diagrams (CMDs) for the measured objects in both fields are 
shown in Figure [3l along with the CMD for NGC 3311 (shown for comparison). In these, 
the GCSs stand out as pairs of vertical sequences. Areas enclosing these vertical sequences 
were chosen as the areas of interest, to eliminate obvious field contamination. Bright limits 
were also chosen so that a range of three magnitudes was initially enclosed (while enclosing 
most of the GCSs, this also allowed the authors to brighten the brightest limits in steps for 
the later analysis discussed in Section [231) . Data points were selected from an initial region 



The catalog can be found at 



http://www. physics. mcmaster.ca/^harris/Databases. html 
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enclosing the GCS vertical sequences within the limits as follows: 

• 0.5 < {g' - i')o < 1.4 and -11.35 < M^, < -8.35 for NGC 5193, 

• 0.25 < {g' - i')o < 1.25 and -11.4 < M^, < -8.4 for IC 4329 and 

• 0.3 < {g' - i')o < 1.2 and -11.3 < M,, < -8.3 for NGC 3311@. 

We note that these selection regions are not identical in colour range, but note again that 
they were chosen with the aim of enclosing the GCS while minimizing field contamination. 
As also noted above, this is likely to be not because of intrinsic differences between the 
galaxies, but rather because of remaining errors in the photometric zeropoints as the different 
galaxies were observed on different nights; but again, the bulk of our later analysis depends 
only on the relative colours and luminosities. The number of data points within this initial 
selection region is shown in column (7) of Tabled We choose to use Mj/ as the GC luminosity 
parameter (as opposed to Mgi) because the values of Mj/ most closely resemble the bolometric 
lumino sity, and because p lotting against the redder colour more accurately portrays the 



MMR (IMieske et al.ll2006l : see also our Section [371]) . The faint-end limits of Mj/>-8.35 and 



-8.4 were estimated visually and adopted to minimize field contamination. The bright-end 
limits of Mj/<-11.4 and -11.35 were only initial and were later increased in steps (see below) 
to see what effect, if any, this would have on the deduced MMRs. 

In Figure HJ we show the radial distribution of the GC candidates around each galaxy 
(see Appendix for a more detailed description of these diagrams, including calculation of the 
GCS specific frequencies). As the CMDs also indicate, the GCS around NGC 5193 is the less 
populous of the two, and is also the more centrally concentrated with a steeper power-law 
fit. 

The first step in the analysis is to test the colour distribution of the GCs for clear 
presence (or absence) of bimodahty. To do this, we used the multimodal stat istical fitting 
package RMIxIfl in the same way as described more fully in lWehner et al.l (120081 ). We carried 



out one-, two- and three- Gaussian fits with the sample divided into colour bins of {g' — i') 
= 0.045, and for comparison we included the previously published NGC 3311 data as well. 
A normal chi-squared test shows that two modes are strongly preferred over only one, but 



^The colour range for NGC 3311 is the same as IWehner et al.l (|2008l ). the fainter limit is also similar, 
and the brighter limit was fixed to form our standard initial enclosure of three magnitudes above the fainter 
limit. 



•^RMIX is publicly available at 



http://www.math.mcmaster.ca/peter/mix/mix.html 
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three modes produced no further improvement with either the hetero- or homoscedastic fits. 
A homoscedastic fit forces the dispersions of both subpopulations to be the same, whereas a 
heteroscedastic fit allows them to be independent. We therefore argue that the two-Gaussian 
fit is the optimal solution, indicating that these galaxies are basically similar to other giant 
ellipticals. The RMIX fits are shown in Figure [5] and all parameters for these fits are shown 
in Table m Columns contain (1) the target's NGC or IC number, (2) the proportion of the 
subpopulation, (3) the error on the proportion, (4) the mean value of the Gaussian fitted 
to the subpopulation (i.e., the peak colour value), (5) the error on the mean, (6) the sigma, 
or distribution value, (7) the error on sigma (not apphcable for the second subpopulation 
because the values of sigma for both subpopulations was forced to be equal), (8) the degrees 
of freedom in the fit, (9) the value, and (10) the reduced value, found by dividing 
by the degrees of freedom. We argue the use of homo- versus heteroscedastic fits for several 
reasons. Significant overlap between the two subpopulations for all three galaxies - NGC 
5193, IC 4329 and NGC 3311 - causes the peaks in the distribution to become obscured. 
Allowing a larger number of free parameters in such a scenario makes an RMIX fit harder 
to converge. When we allow a heteroscedastic fit to the subp o pulat ions of IC 4329, we find 



that the sigma values are equal within errors. IWehner et al.l (120081 ) obtain a similar result 
for NGC 3311. However, for NGC 5193 we find that a heteroscedastic fit has discrepancies 
of a{g' — i') = 0.16±0.01 and 0.10±0.2 for the blue and red subpopulations, respectively. 
With heteroscedastic fits, we also find that the proportion of GCs within each subpopulation 
becomes more unrealistic, especially for NGC 5193: 0.86±0.06 and 0.14±0.06 for the blue 
and red subpopulations in NGC 5193, respectively, 0.61±0.14 and 0.39±0.14 in IC 4329, 
and 0.78±0.04 and 0.22±0.04 in NGC 3311. As we were concerned that these above effects 
were amplified by NGC 5193 GCS's low numbers, we forced all the fits to be homoscedastic. 
Even though the reduced chi-squared value for the two-Gaussian fit on the NGC 3311 GCS 
is the lowest of the different fits tried, we note that it is still unusually high. Various smaller 
bin sizes were tried in the two- Gaussian fit, but the proportions and peak values of the 
subpopulations remained the same. 

The peak values of the red and blue subpopulations for the three Gemini galaxies (in- 
cluding comparison with NGC 3311) are not the same, again because of residual zeropoint 
errors as noted above. However, it is important to note that the differences between the 
colours of the blue and red peaks A{g' — i')(red-biue), are the same among all three galaxies, 
indicating that their metallicity differences are similar. 
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2.3. Analysis of the Mass-Metallicity Relation 



Quite simply, a mass-metallicity relation for either the blue or red GC sequences will 
exist if the mean colour changes systematically with luminosity. The two techniques that 
have previously been used to find it are either to assum e a linear relation between mean 



colour and absolute magri i tude and solye for the slope (IStrader et al.l l2006l : ISpitler et al. 



20061 : iDeGraaff et al.ll2007l : IWehner et al.ll2008l ): or, t o bin the data in magnitude steps and 



epi 

D 



2006 



find the mean colour of each sequenc e in each bin (iHarris et al.l l2006l : iMieske et al 
Wehner et al.ll2008l : [Harris et al.ll2009l ). In this paper, we use the former as the main method 
and compare to the latter. 

The data set for each GCS was sorted into 0.3 magnitude bins in Mj/ so that enough 
bins along the sequence were generated while also maintaining enough GCs per bin to be 
statistically significant. Then to reduce field-star cont amination furt her, we applied an 
approximate representation of the Chauvenet criterion (jParrattI Il96ll ) to each bin. The 
criterion's purpose is to reduce the biasing effect of outlying points on values for the mean 
and standard deviation. We found the mean colour and its uncertainty amean, for each bin 
and discarded outlying points if their colours fell farther from the bin mean than zamean , 
where 



0.85(/o^ion) + 1.12 



and n is the number of points in the bin. 

To get colour-magnitude slopes (i.e., A{g' — i') / AMi/) , we first split the data points into 
red and blue sides with a line at one colour value, although the value of this initial dividing 
line varies from GCS to GCS as discussed below. Figure [6] shows the published data in NGC 
4696, which we use as an example for illustrative purposes. The middle dashed line in this 
figure represents the initial dividing line. The value of the initial di viding line for each GCS 
is chosen to be at the "dip" colour value (IPeng et al.ll2006bl . 



i.e., the colour where a 
note that although the 



GC has equal probability of being red or blue). iPeng et al.l (1200 
mean colours of subpopulations may vary with galaxy luminosity, the dip colour is nearly 
the same. As our data sets have yet to be well-calibrated we can see that our dip colours 
are not all the same: IC 4329 and NGC 3311 have very similar dip colours of (g-i)Q=0.78 
and 0.80, respectively, whereas NGC 5193's dip colour is (g-i)Q=1.02. However, knowing 
that the dip colour values should be constant, we choose the dip colour of the GCS as the 
initial dividing line. We then binned the two data subsets in 0.3 magnitude bins, finding a 
mean for each bin. Finally, we calculated a weighted least-squares (WLSQ) linear fit of bin 
magnitude versus bin colour for each subset (the upper and lower dashed lines in Figure [H]). 
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A line (the middle solid line in Figure [6]) exactly half-way between these two initial 
fits was used to split the Chauvenet-excluded data set a second time, and the process was 
repeated. This second splitting was to compensate for any inequality in the proportions of 
the two subpopulations (see Section 1X21) . The upper and lower solid lines in Figure [H] show 
the final WLSQ fits. Figure [7] shows the CMDs and MMRs for NGC 5193 and IC 4329 
(with only 11 and 12 bins, respectively), where we also added the one-sigma error bars of 
the WLSQ fit to the adopted linear solutions. If we were to eliminate the brightest bin in 
the red slope of IC 4329, we can see that within errors the fitted slope would immediately 
drop to zero. This outlying bin is excluded on similar grounds to the Chauvenet criteria, as 
was a similarly outlying point excluded in NGC 5322. 

As noted above, the upper, brighter limit of the selection region was incrementally 
increased in 0.3 magnitude steps to explore what effect - if any - it had on the colour- 
magnitude slope. We used three methods to do this. For each subpopulation in each GCS 
in our first method, we started with a minimum of ten bins for which we calculated a slope. 
Bins at brighter magnitudes were added until there were no more GCs to bin. With each 
new bin added a new slope was calculated. This gave us up to five slopes per subpopulation 
per GCS. NGC 5193 had so few points that only the slopes for 10 bins and 11 bins could be 
measured (i.e., only one bin was added to the bright end). The two blue A{B — I)q/ AMj^ 
slopes were the same within errors (-0.05±0.03 and -0.05±0.02), as were the red (-0.01±0.05 
and -0.03±0.04). For IC 4329, slopes for 10-13 bins were calculated (i.e., three bins were 
successively added to the bright end). The IC 4329 blue slopes varied between A{B — 
/)o/AM/=-0.11±0.02 and -0.05±0.01, and we found that the red varied between A{B - 
/)o/AM7-=-0.07±0.03 and 0.05±0.02. The values of all blue and red slopes are shown in 
Table [31 The IC 4329 slopes for 12 and 13 bins are shown with and without the outlying 
point. 

The second method used a constant number of ten bins. As a brighter bin was added, 
the faintest bin was dropped and a new slope was calculated. Again, this gave up to five 
slopes per subpopulation per GCS. To look at any possible curvature of the MMR we used 
the IRAF/POLYFIT routine for our third and final method with a polynomial of order two. 
We input both the minimum number of ten bins, and also the maximum number of up to 
fourteen bins in the routine and compared the results. 

With the first two methods, we saw no consistent change on the MMR slopes as brighter 
bins were included. With the final method, most slopes had an insignificant second-order 



'^A{g' — i')/AMii slopes were converted to A{B — I)o/AMj slopes via the relations (jMarastonl 120051 ) in 
Section [3. II to allow comparison to other data sets. Errors are propagated through the equations. 



- 12 - 



component. The remaining few only had a marginal second-order component. We conclude 
that a linear slope for the MMR is a reasonable assumption for our sample of galaxies, and 
that the particular choice of bins (luminosity range) has no major effects on the deduced 
slope. 

To compare all slopes, we used the slope with the maximum number of bins, but exclud- 
ing the faintest and brightest bins. This was done to reduce the effects of field contamination 
(at the faint end) and small numbers (at the bright end). 



3. Comparison and Discussion 



Part of the purpose of our study is to search further for any trends in the slope of the 
MMR with such factors as host galaxy luminosity or environment. To our MMR solutions 
for NGC 5193 and IC 4329, we now add similar res ults from other BC Gs or giant ellipti- 
cals. Previous ly pub lished ones include the eight in [Harris et al.l (120061 ) and NGC 3311 in 



Wehner et al.l (120081). To th ese, we add a further four BCGs - NGCs 7626, 708 a nd 7014, and 



IC 42 96 (IHarris et al.ll2009l ) - observed with HST in {B,I), similar to those in [Harris et al. 
(120061 ). For the total of 15 BCGs listed above. Table [2] lists in successive columns (1) the 



target's NGC or IC number, (2) the galaxy group or cluster, (3) the redshift cz of the galaxy, 
(4) the galaxy luminosity, M^, (5) the f oreground reddenin g, (6) the apparent distance mod- 
ulus (either SBF-based distances from iTonry et al.l (120011 ) were used, as is the case for the 
galaxies marked with (i), or these values were calculated from radial velocities using an 
adopted value of Hq = km s~^ Mpc~^), (7) the number of data points within the initial 
selection region (see Section fl?2\ for details of the selection regions), (8) the total number of 
GCs estimated to be in the system, and (9) the specific frequency, Stv- 



Mieske et al.l (120061 ) conducted an MMR study for 79 early- type galaxies with HST/ ACS 
photometry in ((7475 — 2:850) as part of the Virgo Cluster Survey. Though most of these 
galaxies are dwarfs or ones with small globular cluster populations in which no individual 
MMR solutions can be determined, they combined their sample into four bins grouped by 
galaxy luminosity. In each bin the total GC sample is then large enough for any cumulative 
MMR slope to be found with some confidence. These four bins are compared with the BCG 
data sets. 
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3.1. Comparisons to Other Data 



To put all these studies onto a common system for comparison requires converting 
three different colour indices: ((7475 — zs5q), {B — I), and {g' — i'). In each case we convert 
A{colour) / A{mag) into a scaling exponent p = AlogZ/AlogM. There is currently no estab- 
lished relation between [g' — i') and metallicity for GCs, so we first converted {g' — i') to 
{B — I), although we note that this double conve rsion introduces larger external uncertain- 
ties. With single stellar population models from iMarastoru (120051 ) we calculated the linear 
colour conversion: 



AiB-I) 



(0.53 ±0.12)' 



(2) 



To co nvert to metallicity from (B-I) colours, we used the linear relation as in [Harris et al. 
( l2006l ). which is based on colours and metallicities for Galactic GCs: 



A{B - /)o 
A[Fe/H] 



0.375 ± 0.049 



(3) 



Mieske et al.l (120061 ) obtain their ((7475 — ^sso) colour-metallicity conversion from iPeng et al. 



( l2006bl ). who have a piecewise relation that is based not only on Galactic GCs but also GCs 
in M87 and NGC 4472: 



Aig 



475 



2^850 j 



A[Fe/H] 



0.195 ± 0.025 [0.70 < (c/475 - ^sso) ^ 1-05] 



(4) 



~ ^850 ) 

A[Fe/H] 



0.546 ± 0.069 [1.05 < (^475 - ^sso) < 1-45] 



(5) 



We note that iPeng et al.l (j2006bl ) actually mean log Z rather than [Fe/H]. These two quan- 
tities are different by a factor of [a/Fe], which we assume to be constant so that the slope 
remains the same when converting between colour and metallicity. 

The compilation of results for this entire extended sample of galaxies is shown in Figure 
[HI as a graph of scaling slope p against host galaxy luminosity. Here, we show the MMR 
for each galaxy for both the blue and the red subpopulations. For the iMieske et al.l (120061 ) 
data, each galaxy has two data points plotted. The data set represented by open triangles 
is ((7475 ~ -285o) versus qf475-magnit ude, while the data set represented by solid triangles is 
(5'475 — ^85o) versus 2850-niagnitude. iMieske et al.l argue that the (M^gg^, {g475 — zs5o)) relation 
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is to be preferred, because cross-contamination biasses between the two modes affect the 
(Mg^.g, (5(475 - zs5o)) relation more strongly. We consider only the (M^^^^, (5(475 - zgso)) data 
points (solid triangles) in our discussion. 

For galaxies brighter than My ~ —19, we find that the blue slopes have a mean p— value 
of p = 0.30 ± 0.05. For the red slopes, we find p = —0.1 ± 0.1, which is consistent with zero. 
P- val ues of blue and red slopes are shown in Figure [H The dwarf galaxies in the faintest 



bin of iMieske et al.l (120061 ) have p = 0.11 ± 0.18 for the blue-MMR. However, we note that 
this lower value may be because there are fewer high-luminosity GCs in that bin - even 
after adding all the dwarfs together. Therefore, if the MMR is more prominent at higher 
parent galaxy luminosities, it would affect the dwarf/giant comparison. We find a consistent 
p— value roughly independent of the individual GC luminosity range, although at increasingly 
fainter levels down along the GC sequences the determination becomes less certain because 
of increasing photometric scatter and field contamination. 

Note that our error of Ap=0.05 is an internal error on the slope and does not include 
the transformation errors from {B — I)q to metallicity or from {g' — i') to metallicity, which 
are approximately Ap ^0.02 and 0.2, respectively. These additional external uncertain- 
ties are the inevitable result of combining different photometric systems, but as far as we 
can determine, they do not affect our overall conclusion that p ~ constant for the more 
luminous galaxies. Further, more strictly homogeneous surveys would remove this lingering 
uncertainty. 



3.2. Comparing Different Methods 

We now test whether or not the line splitting method we used above may have actually 
introduced the observed MMR. Suppose we have a situation where the mean colours of both 
red and blue modes stay constant with luminosity, but where the numbers of GCs in the 
two modes might be different. Suppose also (realistically) that the random measurement 
uncertainties increase with magnitude, so that the two modes would overlap more and more 
as we go to fainter magnitudes. If the two subpopulations contained equal numbers of 
GCs, MMRs would not be seen (i.e., they would be consistent with zero) because similar 
numbers of clusters in each mode would contaminate the other mode. However, if the blue 
subpopulation dominated, and the same constant dividing line was used, a positive offset 
for both subpopulations would be seen: more blue clusters would spill into the red side than 
vice versa, so the measured mean colour of both modes would be too blue at the faint end 
of the sample. In this case, we would conclude that the mean colours of both modes became 
redder at brighter magnitudes. Conversely, if the red subpopulation dominated, a negative 
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offset for botli subpopulations would be seen. For tliis reason, we used tlie second iteration 
of tlie line splitting to counter this effect. 

We do not see this contamination-based effect happening in any of the solutions for our 
sample galaxies. For the blue mode, the p— values are consistently positive, whereas for the 
red mode the slopes are slightly negative and consistent with zero. In particular, for the 
red-dominated GCSs (NGCs 1407, 5322, 7049, 3348, 4696ap, 5557 and 7626, and IC 4296) 
the red-mode slopes are all slightly negative or close to zero, which is consistent with the 
above argument. However, in these same GCSs the blue- mode MMRs have an average value 
around p ^ 0.30, which is not consistent with the argument and therefore provides some 
support for the analysis method. 

Another and different factor, which we cannot assess as completely here, is the amount 
of field contamination. Although we selected restricted intervals in colour to minimize the 
contamination, it is notably worse at fainter magnitudes and if the contamination is asym- 
metric i n colour, it could bia s the RMIX fits. An argument against such an effect, however, is 
that the lMieske et al.l (120061 ) Virgo sample is contamination- free (the GCs were identified in- 
dividually) and their mean results are no different from ours at the same galaxy luminosities. 
In addition, the different galaxies in our sample have different degrees of field contamination 
(cf. the individual source papers) but yield very similar MMR slopes. 

We have also compare d our results with na ethods using the KMM bimodal fitting code, 
including six galaxies from [Harris et al.l (120061 ) . and for the three Gemini/GMOS galaxies 
by conducting RMIX fits. (Only six of the eight [Harris et al.l galaxies are used because two 
galaxies only have one bin per subpopulation; each bin contains 200 points.) The Gemini 
data sets were divided into 0.5 magnitude bins, with the brightest two bins in NGC 5193 
and IC 4329 not being included in the fit because of the very low numbers of GCs there. 
Comparisons were made between RMIX and KMM, as shown in Tableland Figure (TUl Six 
of the nine blue-MMRs and a different six of the nine red-MMRs are consistent within errors 
for both methods, further supporting the line splitting method and showing that the codes 
are consistent. 

We find a blue-MMR index at a constant value of 0.30±0.05. However, two galaxies are 
clearly outside that range: NGC 3311, which has a much higher index, and IC 4296 which 
has an index that within the errors is con sistent with z e ro. W e also recall that M49 (NGC 
4472), included in the brightest bin of the lMieske et al.l (120061 ) data set, also does not show 
a blue-MMR. 



- 16 - 



3.3. No Correlation with Specific Frequency 



The fifteen GCSs i n Table [2] can be sp lit into two groups: sparse GCSs, which we define 
in a similar manner to [Harris et al.l (120061 ) to have specific frequency, S'Ar<5; these include 
NGCs 1407, 5322, 3348, 3268, 5557, 7626 and 5193, and ICs 4296 and 4329; and rich GCSs, 
defined to have S'7v>5; these include NGCs 7049, 3258, 4696, 708, 7014 and 3311 (see the 
Appendix in Section [5] for Sn calculations). We find no significant difference in the average 
MMR slope (p— value) between these two groups (see Figure [9]), although we note that NGC 
3311 is a significant outlier. The lack of correlation of p with Sn is consistent with the 
interpretation that the source of the MMR is one more local to the GCs themselves (such 
as self-enrichment) rather than one driven by the large-scale GC formation efficiency. 



4. Summary 

We discuss the photometry for globular cluster systems in fifteen BCGs and use these to 
derive the MMRs for both the red and blue globular cluster subpopulations. Twelve data sets 
were from the HST, and a further three data sets were obtained from Gemini/GMOS with 
{g',i') phot ometry. The e i ght A CS/ WFC data sets a nd on e Gemini data set were previously 



analysed in [Harris et al.l (120061 ) and IWehner et al.l (120081 ) . respectively. For the other two 



CMOS datasets (NGC 5193 and IC 4329) we present the photometry and its results for the 
first time. 

Colour-magnitude slopes were found via a "line splitting" method, and converted to 
an MMR in heavy-element abundance versus mass {Z ~ M^) so that compari sons between 



the H ST and Gemini data, and with the Virgo Cluster Survey galaxy sample (IMieske et al. 



20061 ) data, could be made. We use the composite results to construct a plot of the scal- 
ing exponents - p{blue) and p{red) - versus host galaxy luminosity and also versus specific 
frequency. 

Our work confirms the existence of a blue-MMR with Z^Af°'^°^° °^, independent of host 
galaxy luminosity for any systems brighter than My ^ —19. Individual galaxy-to- galaxy 
differences may still exist (with occasional anomalies such as NGC 4472, whose lack of an 
MMR still remains to be explained) but are still hard to determine unambiguously. For 
the dwarf galaxies fainter than My ~ —19, the blue-MMR slope is consistent with zero, 
though we raise the question that this might be due to the smaller number of GCs at high 
luminosities, making the MMR less noticeable. Whether or not the MMR continues to lower 
GC luminosities with the same slope remains an open question, because of the increasing 
effects of field contamination and photometric measurement scatter. 
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Within the measurement uncertainties, the red-sequence MMR is consistent with a hne 
that includes zero {Z ^ M^^'^^*^'^). We also find no correlation between the p- values of 
slopes and specific frequencies. 

We have tested our line splitting method numerically to confirm that that nature of the 
method does not itself artificially introduce an MMR. In general we find clear agreement 
between the line splitting method and multimodal histogram fitting methods employing 
RMIX and KMM. 



The results of this paper agree partially but not completely with iMieske et al.l (120061 ) . 
who show that for lower host-galaxy luminosities, i.e., dwarf galaxies, the blue-MMR is not 
the same as those for mid-range host-galaxy luminosities. Here, we add many more galaxies 
to the brightest range and find similar relations to those of the mid-range. We note that 
NGCs 3311 and 4296 lie offset from the rest of the galaxies for the blue-MMR in Figure Ei 
Only the dwarf ellipticals stand out as having no clear MMR amongst their globular clusters. 
This suggests that similar GC formation processes occurred in BCGs as occurred in typical 
early-type galaxies. 

The clear existence of a blue-sequence MMR is, in general, consistent with recent models 
in which the higher mean metallicities within more massive GCs are due to self-enrichment 
within the proto-GC: they form in more massi ye potential wells whi ch can hold in a higher 



fraction of supernova-driven enrichment (e.g., iBailin &: Harris! l2009l ) . Within the same in- 



terpretive framework, the lack of a similar trend along the red GC sequence is expected if 
their mean metallicity is due to a higher level of pre-enrichment, above which any later self- 
enrichment would be less important. The independence of p- values versus both host galaxy 
luminosity and specific frequency may favour interpretations where the MMR is due to local 
events in and near where the clusters form, rather than the global environment of the galaxy 
itself. 

We suggest follow-up observations for galaxies as a function of environment, especially 
isolated giant elliptical galaxies. Also, it is statistically difficult to see any MMRs in spiral 
galaxies, but i f they could be stacked together in a similar manner to the dwarf galaxies in 



Mieske et al.l ( 120061 ) then perhaps a trend would become apparent. 
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5. Appendix 

To estimate tlie total number of GCs in tlie Gemini systems, we first counted the number 
of GCs within similar limits of Section 12. 2t 

• 0.5 < {g' - i')o < 1.4 and M^/ < -8.35 for NGC 5193 and 

• 0.25 < {g' - i')o < 1.25 and M^/ < -8.4 for IC 4329. 

These magnitudes were converted to My with the ugriz standard star catalogue by 
E. Wehner referenced in Section 12.11 (i'-V) was plotted against ( g '-i ') and the following 
conversion equation was obtained: 



(i' -V) = -0.69(^' - i') + 0.8238 (6) 

The radial distribution of the objects falling in these magnitude and colour ranges, 
relative to the centre of each galaxy, is shown in Figure HI Here the number density (objects 
per unit area) is plotted against galactocentric distance in linear form, extending out to the 
edges of the GMOS frame. The outer regions of the NGC 5193 field were used to set an 
approximate background level for both fields and this was subtracted from the counts before 
plotting. Then, to estimate the total GC population within the field, we fit a simple power 
law of the form a = (3r~'^ and integrated out to the edge of the image. The power laws, after 
being amplitude adjusted, were a = 313.2 ± 50.4r-2-66±o.22 ^^^^ a = 1.12± i.52r-i-34±o-io 
for NGC 5193 and IC 4329, respectively. As this estimate does not cover the entire radial 
extent of the GCS, we corrected for the field of view in our images. 

We then assumed that for a gE the number of GCs per unit magnitude is describ ed by a 



Gaussian with a mean value of My = -7.33 and a standard deviation of 1.4 magnitudes ([Harris 



20011 ). This allowed us to find the fraction of GCs brighter than the limiting magnitude. By 
dividing the number of GCs above the limiting magnitude by this fraction, the estimate for 
the total number of GCs in the system was finally obtained (i.e., the values in column (8) 
in Table [2]). The specific frequency. Sat is then calculated from 

= Agc10°-'(*'^+^^) (7) 



(IHarris fc van den Berghl Il98ll ). For the remaining galaxies, we either adopted values of 
Sn from the literature or, if there were no previous values, we calculated Sn in a method 
similar to the above without the FOV correction. For HST data, limiting magnitudes were 



converted to My with UBVRI standard star data from iLandoltl (119921 ) . (V-I) was plotted 
against (B-I) and the following conversion equation was obtained: 



{V - I) =0.53{B - I). (8) 
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Table 1. Summary of MMRs in Previously Observed Galaxies 



Galaxy sample 


P 


Brightest GC Luminosity^ 


Reference 


(1) 


(2) 


(3) 


(4) 


8 BCGs 


0.55 


Mj>-13.0 


Harris et al. ("2006") 


NGC 3311 


0.6 


Mj>-13.5 


Wehner et al. (2008) 


Virgo galaxies'' 


0.48±0.08 


Mz>-12.0 


Mieske et al. f2006~) 


M87 


0.48 


M2>-12.2 


Strader et al. (2006) 


NGC 4649 


0.43 


M2>-12.0 


Strader et al. (2006) 


NGC 4472 (M49) 





M2>-11.8 


Strader et al. (2006) 


NGC 5866 


0.11±0.09'= 


Mfl>-10.2 


Cantiello et al. (2007) 


NGC 1533 





M/(8i4)>-10.2 


DeGraaff et al. (2007) 


NGC 4594 (Sombrero Galaxy) 


0.27 


Mfl>-11.6 


Spitlcr ct al. (2006) 



^ Values of cz and extinction from NED were used with an adopted value of 

Ho = 70 km s-l Mpc-l to calculate absolute magnitudes for M87, NGC 4649, NGC 4472, 

NGC 5866, NGC 1533 and NGC 4594. 



^ ^ieske_et_al^ (2006) study 79 early type galaxies which they split into four host-galaxy 
magnitude bins. The MMR cited above is for their brightest bin where the host-galaxy magnitude 
range is -21.7<Ms<-21. 

^antiello_et_al^ (2007) caution that this result is tentative because of the small number of 
GCs (109) used to find these slopes. 
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Table 2. Data for Individual BCG Galaxies 



NGC or IC 


Cluster or Group 


Redshift cz 


E(B-I) 


(m-M)i 


Number of 


Ngc 


Sjv 






(km s-i) 






data points 






(1) 


(2) 


(3) (4) 


(5)- 




(7)- 


(8) 


(9) 


(B, I) photometry from ACS/WFC 


NGC 1407 


Eridanus 


1627' -22.4 


0.16 


32.0 


1046 


2641±303'' 


4.0±1.3'' 


NGC 5322 


CfA 122 


1916' -22.0 


0.03 


32.2 


395 


1600 


2.5 


NGC 7049 


N7049 


1977' -21.8 


0.12 


32.4 


682 


2700 


5.8 


NGC 3348 


CfA 69 


2902=^ -22.2 


0.17 


33.2 


831 


3100 


4.4 


NGC 3258 


Antlia 


3129' -22.1 


0.20 


33.4 


1812 


6000±150<= 


6±2.5° 


NGC 3268 


Antlia 


3084' -22.1 


0.24 


33.4 


1593 


4750±150'= 


3±2« 


NGC 4696 


Cen 30 


3248' -23.0 


0.23 


33.5 


3099 


4100±200f 


6f 


NGC 5557 


CfA 141 


3389'' -22.4 


0.01 


33.4 


841 


3500 


4.1 


NGC 7626 


Pegasus (I) 


3154j -22.4 


0.17 


33.4 


1835 


2833±300S 


3.9s 


B photometry from ACS/WFC, / from WFPC2 


NGC 708 


Abell 262 


4601=* -22.1 


0.21 


34.3 


1282 


4800 


7.2 


NGC 7014 


Abell 3742 


4676=^ -21.8 


0.08 


34.2 


962 


3800 


7.4 


IC 4296 


Abell 3565 


4009=^ -23.4 


0.15 


33.8 


885 


36000 


1.7 


(g, i) photometry from CMOS 


NGC 5193 


Abell 3565 


3988=^ -22.4 


0.11 


33.9 


596 


630±330 


0.70±0.36 


IC 4329 


Abell 3574 


4808=^ -23.1 


0.12 


34.3 


1528 


2970±530 


1.56±0.28 


NGC 3311 


Abell 1060 


3937=" -22.4 


0.15 


34.0 


4936 


16500±2000'' 


12.5±1.5'' 



^ Values were taken from NED. 

^ Values calculated using an adopted value of Hq = 70 km s~^ Mpc 

Points in the initial selection regions as detailed in Scction l2.2l 
Values from 



Perrett et al. (1997) . 
" Bassino et al. f2008') . 
^ Lee fc Geisler (1999) . 
g Sikkema et al. f2006l . 
^ Wehner et al. (2008) . 
' Tonrv et al. ('2001) . 
j Ramella et al. (2002). 
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Table 3. Blue and red A(B-I)o/ AM i slopes with corresponding p values for all bins, 
where p is defined in the relation Z = IVP. Errors are shown as the internal relative error of 
the original slopes (i.e., errors due to the colour transformation are not propagated). The 
IC 4329 slopes for 12 and 13 bins are shown with and without the outlying point 

mentioned in Section [2^31 



NGC 


Number 


Blue 


Blue Slope 


P 


Error 


Red 


Red Slope 


P 


Error 


or IC 


of bins 


Slope'' 


Error" 




on p 


Slope'' 


Error'' 




on p 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 5193 


10 


-0.05 


0.03 


0.31 


0.56 


-0.01 


0.05 


0.05 


0.78 


NGC 5193 


11 


-0.05 


0.02 


0.32 


0.45 


-0.03 


0.04 


0.21 


0.77 


IC 4329 


10 


-0.08 


0.01 


0.51 


0.32 


0.02 


0.02 


-0.14 


0.44 


IC 4329 


11 


-0.05 


0.01 


0.35 


0.33 


0.05 


0.02 


-0.33 


0.44 


IC 4329{-outlier) 


12 


-0.11{-0.111) 


0.02(0.003) 


0.74(0.74) 


0.53(0.30) 


-0.07(0.01) 


0.03(0.02) 


0.48(-0.04) 


0.73(0.36) 


IC 4329{-outlier) 


13 


-0.08{-0.082) 


0.01(0.003) 


0.55(0.55) 


0.42(0.23) 


-0.06(0.01) 


0.03(0.02) 


0.42(-0.06) 


0.64(0.38) 


NGC 3311 


10 


-0.106 


0.003 


0.71 


0.28 


0.006 


0.005 


-0.038 


0.090 


NGC 3311 


11 


-0.103 


0.003 


0.69 


0.27 


0.007 


0.004 


-0.050 


0.080 


NGC 3311 


12 


-0.118 


0.006 


0.79 


0.36 


0.005 


0.003 


-0.033 


0.059 


Values converted from (g-i)o 


to (B-I)o 

















Table 4. RMIX Multimodal Fitting Parameters for Gemini/GMOS Galaxies 



Target 


Pi 


Error 


Mean 


Error 


Sigma 


Error 


DOE 




Reduced 


pr(> X^) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 5193 


0.70 
0.30 


0.05 
0.05 


0.81 
1.11 


0.01 
0.02 


0.144 
0.144 


0.007 

NA 


19 


46.1 


2.4 


2.2x10"^'^ 


IC 4329 


0.68 
0.32 


0.03 
0.03 


0.592 
0.90 


0.009 
0.01 


0.150 
0.150 


0.005 

NA 


18 


69.8 


3.9 


4.8-4 


NGC 3311 


0.57 
0.43 


0.02 
0.02 


0.610 
0.898 


0.005 
0.006 


0.142 
0.142 


0.003 
NA 


19 


221 


11.6 


4.8x10"* 
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Table 5. Comparison of Line Split A{B — I)o/AMj Slopes For Blue and Red GC 
Sequences. Pairs that do not agree within the errors are highhghted in bold. The slopes for 
the line split bins are from the WLSQ fits that have the maximum number of bins minus 

the faintest and brightest bins. 



Galaxy 


Method 


Blue 


Error 


Red 


Error 






Slope 




Slope 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC 1407 


Line split bins 


-0.041 


0.011 


0.0168 


0.0077 




KMM peaks 


0.00 


0.03 


0.03 


0.03 


NOG 3348 


Line split bins 


-0.0517 


0.0081 


0.0148 


0.0065 




KMM peaks 


-0.04 


0.02 


0.027 


0.008 


NGG 3258 


Line split bins 


-0.0495 


0.0027 


0.0083 


0.0064 




KMM peaks 


-0.04 


0.01 


0.00 


0.01 


NGG 3268 


Line split bins 


-0.0571 


0.00.30 


0.0113 


0.0075 




KMM peaks 


-0.061 


0.005 


-0.01 


0.01 


NGG 4696 


Line split bins 


-0.0314 


0.0040 


0.0486 


0.070 


(aperture) 


KMM peaks 


-0.04 


0.01 


-0.04 


0.02 


NGC 5557 


Line split bins 


-0.0498 


0.0054 


-0.0035 


0.0077 




KMM peaks 


-0.030 


0.003 


-0.015 


0.007 


NGG SIQS'^ 


Line split bins 


-0.025 


0.018 


-0.012 


0.032 




RMIX peaks 


-0.11 


0.03 


0.03 


0.09 


IG 4329*^ 


Line split bins 


-0.0400 


0.0089 


-0.033 


0.021 




RMIX peaks 


-0.03 


0.01 


0.03 


0.01 


NGG 3311'' 


Line split bins 


-0.0545 


0.0012 


0.0038 


0.0023 




RMIX peaks 


-0.060 


0.005 


-0.012 


0.007 



Values remain in {g' — i')o/Mj/ 




Fig. 1.— NGC 5193 and IC 4329 Gemini CMOS ^'-band images that have a 5.5' x 5.5' field 
of view and an image scale of 0.146" per pixel, and are from the data set GS-2006A-Q-24. 




Fig. 2. — The NGC 5193 g'-hand and IC 4329 g'-hand data sets are shown with plots of 
instrumental magnitude versus error, and the false identifications within the dashed boxes. 
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g'—i' (instrumental magnitude) g'—i' (instrumental magnitude) 




g I I I I I I I I I I \ I I I I I I I I I I I I I I 

-0.5 0.5 1 1.5 2 
g'—i' (instrumental magnitude) 



Fig. 3. — Colour-magnitude diagrams in {i', g' — i') for the two galaxies imaged in our Gemini 



CMO S program, NGC 5193 and IC 4329, as well as that of NGC 3311 from IWehner et al. 



(120081 ) for comparison. GC candidates are located within the dashed-line box. As discussed 
in Sections 12.21 and 12.31 the brightest limit of the box was incrementally increased in 0.3 
magnitude steps. 
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Fig. 4. — Number of GC candidates per unit area on the sky as a function of projected 
galactocentric radius, after background subtraction. The radial bins used were 50 pixels 
wide, corresponding to 7.3 arcseconds. NGC 5193 and IC 4329 have scales of 3.64 and 2.85 
arcseconds/kpc, respectively. The error bars equal (n°'^/area) for n candidates in each bin. 
The upper horizontal dashed line in each panel shows the adopted background level 0"f,, 
chosen from the outer regions of the NGC 5193 field where the observed number density is 
nearly level. The power-law fits described in the Appendix were determined from the data 
points within the two vertical dashed lines: the inner one represents the point within which 
the detections of the GC candidates become hampered by the increasing fight of the central 
bulge, whereas the outer line represents the largest radius at which a complete annulus fits 
inside the CMOS field. 
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Fig. 5.— The two-Gaussian solutions for the RMIX fits of NGC 5193 and IC 4329 data sets. 
The parameters of these fits are shown in Table |H We note that although the distributions 
look like they could be fit with a skewed unimodal solution, we assume bimodality from 



previ ous studies (e.g.. [Harris et al.ll2006l : iMieske et al.ll2006l : ISpitler et al.ll2006l : IStrader et al. 

2nn6h . 
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Fig. 6. — A CMD for NGC 4696, with an example of initial constant dividing line with a 
value of (B-I)o=1.8 (the middle dashed line) and then a subsequent sloping dividing line (the 
middle solid line). The sloping dividing line is placed exactly half-way between the initial 
WLSQ fits (the upper and lower dashed hues). See Section [2l3l 
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' J : I ■ u i ,| . I , I — m — I — I — I — T — I — I — I — I — I — 
.■; NGC 5193 " 




-8 -g -10 -11 -12 

Fig. 7. — WLSQ fits and their one-sigma errors for the blue and red subpopulations for 
NGC 5193 and IC 4329. The dashed hne shows the initial splitting line and the solid line 
represents the second splitting line. This solid splitting line results in the blue (red) WLSQ 
fit and error lines above (below) it. 
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Fig. 8. — Mean of blue and red MMR indices plotted against host-galaxy luminosity. For our 
data sets (filled circles), the one-sigma errors that were calculated for the colour-magnitude 
slopes are shown as the same relative error on p here. Transformation errors are not included 
here but are approximately Ap ^0.02 and 0.2 for {B — I)o to metallicity transformatio r is and 
{g' — i') to metallicity transformations, respectively. The data points of iMieske et al.l (120061 ) 
are shown here for comparison; {g — z, z) and {g — z, g) as filled and open triangles, respec- 
tively. The errors for their points come from a random resampling around the KMM fits (see 
their Table 1). Note that although we tested 15 GCSs, there are actually 16 points plotted 
above. The two data points at My =-23 are both for NGC 4696, with t he upper (lower) da ta 



point corresponding to the method using PSF (aperture) photometry (IHarris et al.ll2006l ). 
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Fig. 9. — Mean of blue and red MMR indices plotted against the specific frequency, Sn. 
Values and errors for the indices are the same as those in Figure [SI As with Figure M we note 
that although we tested 15 GCSs, there are actually 16 points plot ted above due to both 
PSF and aperture photometry results being included for NGC 4696 (IHarris et al.ll2006l ). 
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-0.15 
-23.5 



-21.5 



Fig. 10. — Comparing line split colour slopes with colour slopes from peaks of KMM and 
RMIX fits as in Table [51 Each GCS has two data points associated with it; hne split slopes 
are on the left (open circles), with KMM/RMIX slopes offset from the true host-galaxy 
magnitude to the right (open squares). 



